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The role of bulk energy in nuclear multifragmentation. 
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Because of thermal expansion and residual interactions, hot nuclear fragments produced in multi- 
fragmentation reactions may have lower nucleon density than the equilibrium density of cold nuclei. 
In terms of liquid-drop model this effect can be taken into account by reducing the bulk energy 
of fragments. We study the influence of this change on fragment yields and isotope distributions 
within the framework of the statistical multifragmentation model. Similarities and differences with 
previously discussed modifications of symmetry and surface energies of nuclei are analyzed. 

PACS numbers: 25.70.Pq , 25.70.Mn , 21.65.+f 



I. INTRODUCTION 



Multifragmentation has been observed in nearly all 
types of high energy nuclear interactions induced by 
hadrons, photons, and heavy ions (see a review 
This is an universal phenomenon occurring when a large 
amount of energy is deposited in a nucleus, and a hot blob 
of nuclear matter is formed. At low excitation energies 
the nuclear system can be treated as a compound nucleus 
, which decays later via evaporation of light particles or 
fission. However, at high excitation energy, or possibly, 
compression during the initial dynamical stage of the re- 
action, this matter blob will expand to the sub-saturation 
densities, where it becomes unstable and breaks up into 
many fragments. As is well known (see e.g. Refs. 
multifragmentation is a fast process, with a character- 
istic time around 100 fm/c. Nevertheless, as shown by 
numerous analysis of experimental data, a high degree of 
equilibration can be reached in these reactions, and sta- 
tistical models are very suitable for description of frag- 
ment yields [!, B 0, 0, B ESEI1 ■ We believe that taking 
multifragmentation into account is crucially important 
for correct description of fragment production in high 
energy reactions. On the other hand multifragmentation 
opens a unique possibility for investigating the phase di- 
agram of nuclear matter at temperatures T w 3 — 8 MeV 
and densities around p w 0.1 — 0.3po (po ~ 0.15 fm~ 3 is 
the normal nuclear density). These conditions are typi- 
cal for the liquid-gas coexistence region. It is interesting 
that similar conditions are realized in stellar matter dur- 
ing the supernova explosions [HI, El| ■ 

In the course of nuclear disintegration hot primary 
fragments are first formed in close vicinity to each other, 
and, therefore, they are still subject to Coulomb and, 
possibly, residual nuclear interactions. It is commonly 
accepted that the liquid-drop description of individual 
nuclei is very successful in nuclear physics. However, in 
a multi-fragment system in the freeze-out volume the pa- 
rameters of the liquid-drop model may change as com- 
pared with those for isolated nuclei. An obvious example 



is the reduction of the fragment Coulomb energy due to 
the presence of other fragments. This effect can be rea- 
sonably evaluated within the Wigner-Seitz approxima- 
tion [If. The Coulomb interaction between the fragments 
can also influence the proton and neutron distributions 
in hot fragments [lj|. Moreover, nuclear interactions pa- 
rameterized as bulk, surface and symmetry energy terms 
in the liquid-drop description of nuclei may change too. 
Possible modifications of surface and symmetry energies 
of primary fragments, and constraints from relevant ex- 
perimental data, were analyzed in the previous works 
[H, EE E3, El, EH- :t was found that the symmetry 
energy of hot fragments drops significantly in the freeze- 
out volume, and the surface energy can be considerably 
modified at high temperatures. In this paper we investi- 
gate possible changes of the nuclear bulk energy in pri- 
mary fragments, and how this may affect their produc- 
tion in multifragmentation reactions. We will show that 
these effects may be quite essential for explaining some 
key observables, and they should be included in realistic 
models. 

2. Statistical description of nuclear multifrag- 
mentation 

All dynamical models used for description of the initial 
stage of the reaction lead to the conclusion that after a 
time interval of few tens fm/c, when fast particles leave 
the system, the evolution of the remaining nuclear system 
changes its character. Because of intensive interactions 
between nucleons the system evolves toward statistical 
equilibrium. At later times the hot nuclear residue ex- 
pands and breaks-up into hot primary fragments. The 
Statistical Multifragmentation Model (SMM) is based on 
the assumption of statistical equilibrium between pro- 
duced fragments in a low-density freeze-out volume [l|. 
We believe that at this point the chemical equilibrium 
is established, i.e., the baryon composition (mass and 
charge) of primary fragments is fixed. However, the frag- 
ments can still interact with other nuclear species via the 
Coulomb and nuclear mean fields. Hence their energies 
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and densities may be affected by these residual interac- 
tions. All breakup channels composed of nucleons and 
excited fragments are considered, and the conservation 
of mass, charge, momentum and energy is taken into ac- 
count. An advantage of the model is that the formation 
of a compound nucleus is included as one of the channels. 
This allows for a smooth transition from the decay via 
evaporation and fission at low excitation energies [2| to 
the multifragmentation at high excitations. In the micro- 
canonical treatment [3, [2(| the statistical weight of the 
decay channel j is given by W- } oc exp Sj , where Sj is the 
entropy of the system in channel j which is a function 
of the excitation energy E x , mass number Aq, charge Zq 
and other global parameters of the source. After for- 
mation in the freeze-out volume, the fragments propa- 
gate independently in their mutual Coulomb field and 
undergo secondary decays. The deexcitation of the hot 
primary fragments proceeds via evaporation, fission, or 
Fermi-breakup [21] . 

In the SMM light fragments with mass number A < 4 
and charge Z < 2 are considered as structure-less parti- 
cles (nuclear gas) with masses and spins taken from the 
nuclear tables. Only translational degrees of freedom of 
these particles contribute to the entropy of the system. 
Fragments with A > 4 are treated as heated nuclear liq- 
uid drops, and their individual free energies Faz are pa- 
rameterized as a sum of the bulk, surface, Coulomb and 
symmetry energy terms: 
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In this standard expression F^ z = (—Wo — T 2 /eo)A is 
the bulk energy term including contribution of internal 
excitations controlled by the level-density parameter eo, 
and Wo = 16 MeV is the binding energy of infinite nu- 
clear matter. F% z = BoA 2 / 3 ((T 2 - T 2 )/(T 2 + T 2 )) 5 / 4 is 
the surface energy term, where Bo = 18 MeV is the sur- 
face coefficient at T = 0, and T c = 18 MeV is the critical 
temperature of infinite nuclear matter. The Coulomb 
energy is E AZ = cZ 2 /A 1 ^ 3 , where c is the Coulomb 
parameter obtained in the Wigner-Seitz approximation, 
c = (3/5)(e 2 /r )(l — (p/po) 1 ^ 3 ), where e is the pro- 
ton charge, ro=1.17 fm, and the last factor describes 
the screening effect due to presence of other fragments. 
E A v z l — j(A — 2Z) 2 /A is the symmetry energy term, 
where 7 = 25 MeV is the symmetry energy coefficient. 
These parameters are taken from Bethe-Weizsacker for- 
mula and correspond to the isolated fragments with nor- 
mal nuclear density. This assumption has been proven 
to be quite successful in many applications. However, a 
realistic treatment of primary fragments in the freeze-out 
volume may require certain modifications of the liquid- 
drop parameters as indicated by experimental data. 

In the grand canonical treatment of the SMM [22j , af- 
ter integrating out translational degrees of freedom, one 
can write the mean multiplicity of nuclear fragments with 
A and Z as 
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Here gAz is the ground-state degeneracy factor of species 

1/2 

(A,Z), At = (2ttH 2 /watT) is the nucleon thermal 
wavelength, and m^r ~ 939 MeV is the average nucleon 
mass. Vf is the free volume available for the translational 
motion of fragments. The chemical potentials p and v are 
found from the mass and charge constraints: 
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As was demonstrated by numerous comparisons of the 
SMM with various experiments, the model describes data 
very well (see, e.g., Refs. @, @, 0, 1, !, M, EH ) ■ This con- 
firms that the statistical approach with liquid-drop de- 
scription of individual fragments provides adequate treat- 
ment of the multifragmentation process. This also justi- 
fies the application of the statistical approach for investi- 
gating the liquid-gas phase transition in nuclear systems 

Mm 

3. Modifications of properties of primary frag- 
ments 

In recent years several new analyzes of experimen- 
tal data with statistical models, related to the nuclear 
isospin in multi frag mentation reactions have been per- 
formed [H, [HI \M [H[. They conclude that modifica- 
tions of the liquid-drop parameters of hot fragments pro- 
duced in the freeze-out volume are needed to explain the 
data. It was suggested that this can happen because of 
a new physical environment where fragments are formed, 
in particular, since they are surrounded by nucleons and 
other hot fragments. The residual interactions may lead 
to energy and density changes which can effectively be ex- 
plained by a modification of the macroscopic nuclear pa- 
rameters. The symmetry energy coefficient 7 was inves- 
tigated in several independent experiments H^, [H, [n| , 
which used both the isoscaling phenomenon [251 ] and iso- 
tope distributions of fragments. It is important that all 
experiments come to the conclusion that the coefficient 
7 drops from about 25 MeV, known for isolated cold 
nuclei, down to ~ 15 MeV for hot primary fragments 
at multifragmentation. The same results were extracted 
also from analysis of mean neutron content of fragments 
[T3L |]~9l| . One of the aims of future experiments is to verify 
this conclusion. 

A recent analysis of the ALADIN data has revealed 
modifications in the nuclear surface properties too [l5| . 
There the neutron-to-proton (N/Z) dependence of the 
surface energy was analyzed for different event classes 
corresponding to different excitation energies. At low 
excitation energies, corresponding to the onset of multi- 
fragmentation, the surface energy follows the trend pre- 
dicted by the standard liquid-drop model, i.e., it de- 
creases with N/Z. This trend is usually explained by the 
contribution of the surface part of the symmetry energy. 
In the region of developed multifragmentation (tempera- 
tures T ss5-6 MeV), where intermediate-mass fragments 
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(IMF: Z = 3 — 20) are mostly produced, the surface en- 
ergy coefficient becomes nearly independent of the N/Z 
ratio. Taking into account this result we conclude that 
subdivision of the total symmetry energy into volume and 
surface parts becomes irrelevant at multifragmentation 
conditions. We point out also that similar conclusions 
about changing surface properties of nuclei were obtained 
within the dynamical AMD model [26|. Therefore, the 
observed decrease of the symmetry energy of fragments 
should not be related to the increase of the total surface 
of fragments at multifragmentation. It is more likely that 
this is an indication of the medium modification of 7 co- 
efficient in the system of many fragments. 

We should point out that there are attempts to ex- 
plain the isospin related observables within the mean- 
field models, as a result of reduced mean density of the 
excited nuclear system [27l [28| ■ As is well known, the cor- 
responding equations of state predict decreasing symme- 
try energy at subnuclear densities. However, within this 
approach the fragment formation is not considered, and 
it is assumed that the produced fragments keep a " mem- 
ory" about symmetry energy of a low-density system. In 
our opinion this assumption is not very convincing, since 
the fragment formation process (e.g., associated with the 
liquid-gas phase transition) may drastically change prop- 
erties of the system, and especially properties of primary 
fragments. We believe that it is more reasonable to as- 
sume that primary fragments produced in an expanded 
freeze-out volume have a reduced density, and this can 
lead to reduction of their symmetry energy observed in 
the experiments fvA fl8l. Ii~9j . 

Until now the nucleon density of individual fragments 
was not explicitly considered in the SMM calculations, 
and, on default, it was assumed that all fragments have 
normal nuclear density po- However, as a result of resid- 
ual interaction between the fragments and their thermal 
expansion the density of fragments (pf) may in fact be- 
come smaller than po- This will lead to a reduction of 
the bulk binding energy Wo- For estimation of this effect 
we use the following expression for the bulk energy as a 
function of the nucleon density: 

^(P /)S -W ( Pf ) = -W + ^-Bl^2l, (3) 

where K w 260 MeV is the nuclear compressibility mod- 
ulus. There are natural limits for possible reduction of 
the average density of fragments. Due to the spinodal 
instabilities and Coulomb fluctuations at the subnuclear 
densities p < (0.6 — 0.7)/?o, uniform nuclear matter be- 
comes unstable, and the nuclear "pasta" phases are pro- 
duced [1^, [3^, [H| . By this reason, we take Pf ~ 0.6po as 
the minimum possible nuclear density in individual frag- 
ments. One can see from eq. ([3]), that this corresponds to 
decreasing Wq from 16 MeV at pf = p down to Wq » 14 
MeV. In the following analysis wc allow for even smaller 
Wq, in order to perform a complete investigation of this 
effect. 

At the last stage of the multifragmentation process hot 



primary fragments undergo deexcitation and propagate 
in the mutual Coulomb field. As was demonstrated in 
many works (see, e.g., [l|) this stage is very important 
for correct calculations of final yields of fragments. In the 
beginning of the deexcitation the hot fragments are still 
surrounded by other species, and, therefore, their mod- 
ified properties should be taken into account. As far as 
we know, only one evaporation code was designed, which 
takes into account the modified properties of fragments 
in their de-excitation. It was developed in Refs. [la.[32T|. 
where modifications of symmetry energy were explicitly 
considered. In the present analysis we use the same pre- 
scription. Namely, we start from the modified bulk en- 
ergies of hot nuclei and restore their normal properties 
by the end of the evaporation cascade. In actual calcula- 
tions we have used a simple interpolation between these 
two limiting states. The energy and momentum conser- 
vation laws were fulfilled in the course of this evaporation 
process. 

We emphasize that the treatment of this de-excitation 
stage should be consistent with the physical properties 
of primary fragments. For example, a failure to describe 
the experimental isoscaling data [33| , by using a dynam- 
ical model for primary fragment formation, may be re- 
lated to the fact that the sequential evaporation from 
primary fragments was included without paying atten- 
tion to their reduced densities and, consequently, to their 
lower symmetry energies. As was demonstrated in Refs. 

, [l8j , this effect can influence essentially the observed 
isotope distributions, since separation energies of neu- 
trons and protons are changed when compared to the iso- 
lated low-excited nuclei. The transition to the secondary 
de-excitation can be consistently controlled in statistical 
description of fragment formation, by using a general- 
ized evaporation prescription [l6l |. Moreover, we think 
that the final conclusions about properties of primary 
fragments can be obtained only after a many-component 
analysis of experimental data. Besides the isotope infor- 
mation and isoscaling observable this should include the 
corresponding fragment charge distributions, IMF mul- 
tiplicities, temperatures, and the other relevant charac- 
teristics (examples of such an approach are presented in 
Refs. dSSEii). 

4. Influence of the bulk energy on multifrag- 
mentation characteristics 

We have used the SMM to simulate multifragmentation 
of the gold source with excitation energies in the range 
of 2-12 MeV per nucleon at different values of the bulk 
energy coefficient Wo- For simplicity, wc fix the freeze- 
out density at l/3po as was done also in many previous 
SMM studies. In Figure 1 we show some characteristics 
of the hot fragments: effective temperatures, mass num- 
ber of the largest fragment, average multiplicity of IMFs 
produced in the freeze-out volume. The effective temper- 
ature T e ff is calculated from the energy balance in the 
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system, in the same way as it was done in Ref. [1 
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E x (MeV/nucleon) 

FIG. 1: Effective temperature, reduced mass number of the 
maximum fragments and average multiplicity of the interme- 
diate mass fragments as a function of the excitation energy 
at different bulk energy coefficients Wo = 12, 14, 16 MeV are 
shown, respectively, in the top, middle and the bottom panels. 



One can see that decreasing the bulk energy does not 
change qualitatively the general picture of multifragmen- 
tation. Nevertheless, some new important features ap- 
pear. For example, the back-bending of the caloric curve 
become more pronounced, i.e., multifragmentation reac- 
tions become more endothermic. Such behavior of T e ff 
can be interpreted as a manifestation of the negative heat 
capacity. An interesting result is that the number of 
IMFs and the mass of large fragments at freeze-out be- 
come smaller, while the number of light particles (Fig. 2) 
is increased considerably. From Fig. 2 one can see a dras- 
tic rise of the yields of a-particles and clusters with A < 4 
at reduced Wo- The reason is that we take into account 
fragments with A < 4 as 'gas' particles with their table 
binding energies, while binding energies of other frag- 
ments become smaller. We believe that a-particles have 
normal properties in the low density freeze-out volume. 
This assumption is consistent with many experimental 
data. For example, by interpretation of fragment yields 
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FIG. 2: Yields of the neutrons, protons, tritons and alpha 
particles in the hot freeze-out volume as a function of the 
excitation energy at different bulk energy coefficients. 



in direct knock-out reactions on nuclei, a- clusters are 
assumed to pre-exist even at the normal nuclear density, 
which is much higher than the freeze-out density at mul- 
tifragmentation |35j. We suggest that this effect may 
explain an excess of a- particle yields over SMM predic- 
tions observed in emulsion data [36j . 

It is important that all new trends in the fragment 
yields caused by changing the bulk energy survive after 
the secondary deexcitation. Figs. 3 and 4 show the fi- 
nal yields of cold fragments. The number of emitted a- 
particles even increases, since their evaporation becomes 
more probable at the first stages of deexcitation, in ad- 
dition to the increased production at the break-up of the 
source (see Fig. 4). This fact leads to the decrease of the 
so-called He-Li temperature (see Fig. 3), which is one 
of the useful experimental observable for the calorimetry 
of multifragmentation processes [13, [H[ ■ On the other 
hand, decreasing Wo causes a drop of neutron and pro- 
ton yields, since they are now bound in the light clusters. 
For the same reason, the maximum IMF multiplicity be- 
comes lower. 

In order to give a detailed picture of nuclear break-up, 



5 




E x (MeV/nucleon) 

FIG. 3: Characteristics of the final cold fragments (after sec- 
ondary deexcitation) produced from Au sources. The nota- 
tions are the same as in Fig. 1, here in the top panel variations 
of the isotopic Helium-Lithium temperatures The-li are also 
shown. 



in Figs. 5 and 6 we present the distributions of produced 
fragments in the full mass range, before and after the 
secondary deexcitation. We choose the range of excita- 
tion energies where the transition from the 'U-shaped' to 
the power-law distributions is taking place. It is obvious 
that decreasing the bulk energy favours multifragmenta- 
tion, since less bound largest clusters are destroyed in 
favor of light clusters and nucleons. However, the yields 
of IMFs in this crucial region change very little with Wo 
(as one can also see from Figs. 1 and 3), and the general 
shape of the IMF distributions does not change. By us- 
ing a well-known power-law parametrization A~ T of IMF 
yields we obtain that the r remains practically the same 
for the considered variations of Wq. This result is quite 
understandable: since the Wq parameter enters binding 
energies of all IMF, their relative yields do not change. 
For this reason, the previously reported results concern- 
ing evolution of the isospin-dependent contribution to the 
surface energy with excitation energy of the system [l5| , 
which are based on consideration of the IMF yields, re- 
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FIG. 4: Yields of the neutrons, protons, tritons and alpha 
particles for final cold fragments, after secondary deexcita- 
tion. 



main valid in the case of reduced bulk energy too. 

Much attention is now paid to the isotope production 
in multifragmentation reactions, because of its connec- 
tion to properties of neutron-rich nuclear matter and to 
astrophysical applications [HI . We have investigated this 
problem, and Fig. 7 demonstrates isotope distributions of 
fragments with Z — 10 before (in the freeze-out volume) 
and after the secondary deexcitation (at infinity). At 
smaller bulk energies the hot fragments become a little 
bit more neutron rich. This is explained by enhanced 
production of the symmetric a clusters, leading to a neu- 
tron enrichment of the remaining nuclear matter. How- 
ever, this trend is more evident for the cold fragments, 
and we have two contributions to this effect. The first 
one is caused by the considerable probability of a parti- 
cles emission in the beginning of the evaporation cascade, 
due to the modified bulk energy of heavier nuclei. The 
second one is related to the fact that the initial tempera- 
ture is lower, and, therefore, the evaporational evolution 
of nuclei towards the /3-stability line is ceased at rela- 
tively large neutron excess. This possibility for obtain- 
ing neutron-rich isotopes should be considered in future 
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FIG. 5: Mass distributions of hot fragments produced from 
Au sources at excitation energies of 4, 5 and 6 MeV per nu- 
cleon for different bulk energy coefficients. 
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FIG. 6: Mass distributions of the cold fragments produced 
from Au sources at excitation energies of 4, 5 and 6 MeV per 
nucleon. 



analyzes of data alongside with possible reduction of the 
symmetry energy of fragments [2, [l7], El, El, HH . 

On the other hand, these two possibilities may com- 
plement each other, since an expansion (i.e., a reduced 
density) of hot fragments should lead to decreasing their 
symmetry energy too [3!| . The bulk symmetry energy as 
a function of density is usually parameterized as 

7 (p / )oc 7 ( /90 )(^r , (4) 
Po 

where the exponent n ranges between 0.5 and 1.5 depend- 
ing on the model assumptions on isospin-dependent nu- 
clear interactions. Recent analyses of experimental data 
have shown that around on the values of n may lie in 
the interval between 0.7 [5o] and 0.9 [4l|. Thus, in or- 
der to explain a considerable drop of 40% in 7, which 
fits experimental observations, we need to decrease pf 
to (0.45 — 0.55)po- As we have noted previously, this 
is already lower than the threshold for appearance of 
different "pasta" phases [29|, [3(| [Mj]. Also, the realis- 
tic Hartree-Fock and Thomas-Fermi calculations for hot 
isolated compound nuclei predict a rather moderate ex- 



pansion at temperatures of 5 - 6 MeV, not more than 
10 - 20% [HI, [43] , depending on the nuclear forces used. 
Therefore, the maximum reduction of the fragment den- 
sity, which can be approximately considered, is around 
30-40%. In this situation one should find additional 
mechanisms to explain the observed modification of the 
symmetry energy. 

5. Conclusions 

In this paper we proceed further with the investiga- 
tion of possible modifications of nuclear properties in sur- 
rounding of other nuclei, which can still interact with nu- 
clear and electromagnetic forces. This subject is relevant 
for nuclear matter at low density, which breaks-up into 
fragments. This state of matter can also be considered as 
a mixed phase of the nuclear liquid-gas phase transition. 
It is expected that this transition takes place in such as- 
trophysical processes as collapses of massive stars and su- 
pernova explosions. In terrestrial laboratories this phys- 
ical phenomenon can experimentally be studied in mul- 
tifragmentation reactions induced by hadrons and heavy 
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FIG. 7: Mass distributions of hot and cold fragments with 
the charge number Z=10 at the excitation energy of 5 MeV 
per nucleon. 



ions of intermediate energies. Actually, the multifrag- 



mentation is a unique tool, which extends our conven- 
tional study of isolated compound nuclei to the nuclei 
embedded in the environment of many nuclear species. 

In recent years interesting new analyzes of multifrag- 
mentation data have appeared which indicate a sig- 
nificant modification of properties of nuclei in such a 
medium. They include: decrease of the symmetry energy 
of nuclei, disappearance of an isospin-dependent contri- 
bution to the surface energy, and possible expansion of 
hot nuclei. In the present work we have investigated the 
influence of decreasing bulk energy on the fragment pro- 
duction. We have found that this effect should lead to 
an enhanced production of a clusters. In the case of a 
neutron-rich system, this should result in larger neutron 
richness of heavy fragments. Another important finding 
is that, despite of disintegration of the system into larger 
number of fragments, the shape of the IMF mass dis- 
tributions does not change with decreasing bulk energy. 
Also we have found a moderate modification of the caloric 
curve, namely a more pronounced back-bending, which 
indicates a negative heat capacity, and a drop of the ap- 
parent helium-lithium temperature. It is important that 
all findings are consistent with previous results on mod- 
ification of fragment properties extracted from experi- 
mental data. Moreover, the expansion of hot fragments 
may provide an additional contribution to some effects 
observed earlier, such as decrease of their symmetry en- 
ergy, and increase of neutron richness of cold fragments. 
We foresee new analyzes of experimental data in order to 
clarify contributions of different mechanisms leading to 
the modification of fragment properties in multi-particle 
environment. 
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